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Testing the effects of opacity and the chemical mixture on the 
excitation of pulsations in B stars of the Magellanic Clouds 
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ABSTRACT 

The B-type pulsators known as /? Cephei and Slowly Pulsating B (SPB) stars present pulsa- 
tions driven by the k mechanism, which operates thanks to an opacity bump due to the iron 
group elements. In low-metallicity environments such as the Magellanic Clouds, /3 Cep and 
SPB pulsations are not expected. Nevertheless, recent observations show evidence for the 
presence of B-type pulsator candidates in both galaxies. 

We seek an explanation for the excitation of yS Cep and SPB modes in those galaxies by ex- 
amining basic input physics in stellar modelling: i) the specific metal mixture of B-type stars 
in the Magellanic Clouds; ii) the role of a potential underestimation of stellar opacities. 
We first derive the present-day chemical mixtures of B-type stars in the Magellanic Clouds. 
Then, we compute stellar models for that metal mixture and perform a non-adiabatic analysis 
of these models. In a second approach, we simulate parametric enhancements of stellar opac- 
ities due to different iron group elements. We then study their effects in models of B stars and 
their stability. 

We find that adopting a representative chemical mixture of B stars in the Small Magellanic 
Cloud cannot explain the presence of B-type pulsators there. An increase of the opacity in the 
region of the iron-group bump could drive B-type pulsations, but only if this increase occurs 
at the temperature corresponding to the maximum contribution of Ni to this opacity bump. 
We recommend an accurate computation of Ni opacity to understand B-type pulsators in the 
Small Magellanic Cloud, as well as the frequency domain observed in some Galactic hybrid 
/3 Cep-SPB stars. 
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1 INTRODUCTION 

j8 Cephei and SPB pulsators are main sequence (hereafter MS) B 
stars presenting self-driven oscillations. The former present low- 
order acoustic and gravity modes (respectively p and g modes), 
while the latter show high-order g modes. yS Cep stars have masses 
from 7 to 20 Mp, with typical periods between 2 and 8 hours 
JStankov &"Han dler 2005). SPB masses r ange from 2.5 to 8 Mo, 
with typical periods between 1 and 3 days ( I Waelkensll 1 99 lb . These 
oscillations are excited through the k mechanism, activated by the 
presence of an opacit y bump at T ^ 200, 000 K due t o the iron 
group elements (see Moskalik & Dziembowski 1992|; ICox et al.l 
ll992l :lDzie mbowski. Mos kalik & Pamvatnvkh 1993). 

The opacity profile, shaped by the iron group elements, de- 
termines the limits of the instability strip of B-type pulsators and 
the frequency domain of their oscillations. Individual abundances 
of the iron group elements (Fe, Ni, Cr and Mn) and their dis- 
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tribution inside the star, as well as their contribution to stellar 
opacity are essential to understanding the seismic properties of 
B-type pulsators. Indeed, the nature of yS Cep pulsations could 
only be understood onc e new op acity computations by OPAL 
(Iglesias & Rogers '1996) and OP (Seaton et a"i]|l994h . including 
the contribution of a huge number of bound-bound transitions from 
the iron group elements, were available. These opacity computa- 
tions were then able to account for most of the obser vational data 
for/? Cep and SPB pulsators available at the time (e.g. IPamvatnvkhl 
1 19991) and succes sfully fit asteroseismic observables of j3 Cep stars 
such as 16 Lac JDziembowski & lerzvkiewica 19961: Thoul et al.l 
l2003h. HP 129929 JAerts et alj2003l : lDupret et al.l2004 or j8 CMa 
I Mazumdar et al.ll2006h . 



With progress in observational techniques, more and more fre- 
quencies were observed in /? Cep and SPB stars while objects pre- 
senting both low-order p and high-order g modes (hybrid pulsators) 
were found. As a consequence, new discrepancies between obser- 
vational data and theoretical predictions appeared. Models were 
not able to account for oscillation spectra of hybrid pulsators (e.g. 
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12 Lac and v Eri: Pamvatnvkh. Handler & Dziembowskil |2004|: 



lAusseloos et alj l2004h . that is for the simultaneous detection of 
modes spanning a large frequency domain, including SPB and /? 
Cep pulsations. 

In the last decade, two events, t he update of the solar mixture 
JAsplund, Grevesse & Sauvall I2OO5, hereafter AG S05) and new 
opacity computations by OP JBadnell et alj|2005r) . led to impor- 
tant consequences for the properties of B-type pulsators. The re- 
vision of the solar mixture decreased the abundances of C, N, O 
and Ne by ~ 30% with respect to the "standard"solar mixture 
dOrevesse & Noelsl [19931 hereafter GN93), while the Fe one did 
not change. For a fixed metallicity (Z), the mass fraction of Fe is 
then 25% larger in the new mixture than in GN93. On the other 
hand, the OP team improved opacity computations JBadnell et al.l 
I2OO5I) by including new inner-shell atomic data for 13 elements 
and by complementing previous data with new inter-combination 
lines and photo-ionization cross-sections. In particular, this led to 
an 18% increase of the Rosseland mean opacity in the region of 
the iron opacity bump. As a result, for the structure of typical 
P Cep st ars, the iron opacity bump is modified by up to ~ 30% (see 
Fig. 4 oflBadnell et ahOOOl see also'Mi eUo. Montalban & Dupretl 
l2007ah compared to what would be obtained with the coimnonly 
used OPAL opaci t ies jlglesias & Rogerslll99q) . 

iMiglio et al.l d2007ah studied the eff'ects of the solar mixture 
(GN93 vs AGS05) and the opacity data (OPAL vs OP) on mode 
stability in MS B stars and found that the most favourable case 
for the excitation of /3 Cep and SPB modes was obtained by com- 
bining the AGS05 mixture with the new OP data. The relative in- 
crease, at a given Z, of the mass fraction of iron group elements 
in the AGS05 mixture yields wider instability bands than does 
the GN93 mixture. The hot wing of the opacity bump in OP data 
also extends the blue edge (hotter limit) of the instability strip in 
comparison with OPAL. Finally, the OP opacity profile leads to a 
larger frequency domain of excited modes and to a greater num- 
ber of hybrid pulsators, as SPB modes are more easil y excited 
(see also iMiglio. Mo ntalban & Dupretl l2007bl : |Pamvatnvki]l2007l : 
IZdravkov & Pamvatnvkh.2008i) . 

The combination of the AGS05 mixture and the up- 
dated OP opacities has no w been used in several studies. 
iDziembowski & Pamvatnvkhl |2008 ) obtained a lar ger frequency 
domai n of excited modes in modelling v Eri while iBriguet et al.l 
( I2OO7I) had promising results for the /? Cep star 6 Oph. How- 
ever, the lowest frequency modes detected in v Eri still can- 
not b e theoretically reproduced (IDziembowski & Pamvatnykhl 
120080 . Additional seismic studies show that OPAL is better 
at reproducing the n on-adiabatic properties of (3 Cep modes 
in 6 Oph and v Eri dDaszvnska-Daszkiewicz & WalczakI l2009l 
|201(]|). whereas OP is preferable for an SPB mode in v Eri 
I Daszyiiska-Daszkiewicz & WalczakI I2OI Of) . The excitation of the 
lowest frequenci es of 12 Lac are however still not explained by the- 
oretic al models ( iDziembowski & Pamvatnvklill2008l : JDesmet et al.l 
I2OO9I) . 

The disagreement between observations and theoretical pre- 
dictions based on updated input physics is even stronger when low- 
metallicity environm ents are considered. T heoretical computations 
with AGS05 and OP JMiglio et alj|2007ij) predict SPBs andyS Cep 
stars down to respectively Z = 0.005 and Z = 0.01. Therefore, these 
pulsators should not be present in the Small Magellanic Cloud 
(SMC) and only rarely in the L arge Magellanic Cloud (LMC), 
galaxies for which JBuchlen ( I2OO8I) , from a study of Cepheid prop- 
erties, derived average metallicities of Z = 0.0027 for the SMC and 
Z = 0.0046 for the LMC. Nevertheless. JKolaczkowski et al.l IJOO^ 



announced the discovery of 151 candidates (92 yS Cep stars and 59 
SPBs) in t he LMC, and 17 (6 /3 Cep stars and 1 1 SPBs) in the SMC. 
Moreover, Sarro et al.l (12009 ') found 63 candidates (48 yS Cep stars 
and 15 SPBs) in the LMC and 43 (23 /? Cep stars and 20 SPBs) in 
the SMC. B-type pulsators have also been detected in SMC s tellar 
clusters: 29 SPB candidates in NGC 371 ( JKaroflFet al.ll2008l), and 
1 jgC ep candidate and 8 SPB candidates in NGC 330 JDiago et all 
I2OO8I) . 

All that calls into question some of the input physics used in 
modelling opacity profiles and consequently the stability of pulsa- 
tion modes. It is generally assumed that for population I stars, the 
distribution of element abundances in the stellar matter is the same 
as in the Sun, although scaled to th eir metallicity (Z). However, 
as sho wn by Mi glio et aU (.2007a) and lMontalban. Miglio & Morell 
( 12009!) ■ changing the pattern of element abundances (even with the 
same value of Z) can have important implications for the stability 
analysis of B-type pulsators. Hence, we wonder whether the solar 
mixture is representative of MS B stars in the MCs. 

The problem could also arise from opacity calculations which 
may underestimate contributions in some particular temperature 
and density domain. This possibility has already been invoked to 
solve old and recent disagreements between theore tical and ob- 
servational seismic properties of MS pulsat o rs (e.g. Isimonll 19821: 
JMontalban et al.' I2OO4I; iBahcall et al.l boO^; JMontalban & Mighd 
2008; Dziembowski & Pamyatnvkh 2008), as well a s of evolved B- 
type pulsators (lefferv & Saio 2006). In particular, jjefl^erv & Said 
( 12006) showed that increasing the opacity due to Ni could help with 
fitting the observational subdwarf B (sdB) instability strip. Fortu- 
nately, these remaining problems have motivated different teams 
to undertake experiments trying to reproduce stellar plasma con- 
ditions in a laboratory to measure monochromatic opacities (e.g. 
JLoiselet al.i 2009) . The first comparisons with calculations made 
by OP and other plasma opacity codes, including attempts to im- 
prove theoretical calculations, reveal discrepancies, in particular for 
the Fe and Ni monochromatic opacities (Turck-Chieze et al. 2011; 
Gilles et al. 2011, and references therein), suggesting that there is 
still room for improvement in opacity data, with important conse- 
quences for B-type pulsators. 

Although we are aware that non-standard processes might 
have some effects on B star pulsations, our aim in this paper is to 
address the role of standard input physics on the driving of pulsa- 
tions in MS B stars. We begin with detailing the chemical mixtures 
representative of B stars in the MCs in the following section. In 
the next section, we then consider the pulsation stability of B star 
models computed with the SMC chemical mixture. In Sect. [4] we 
analyse the effects of an ad hoc local change in the opacity, to iden- 
tify the changes that could explain the presence of B-type pulsators 
in the SMC. A discussion and a conclusion end the paper 



2 CHEMICAL MIXTURE OF YOUNG STELLAR 
OBJECTS IN THE SMC 

The detailed chemical mixture and opacity data are two basic fac- 
tors that can affect the driving of pulsations in/? Cep stars and SPBs. 
The Magellanic Clouds are a very convenient place to observe not 
only the effects due to a different chemical mixture but also those 
coming from a low metallicity. The results of Miglio et al. (2007a|) 
predict a limited number of SPBs and no /3 Cep stars in the MCs. 
These results are challenged by recent detection of a significant 
number of B-type pulsator candidates in the MCs, and particularly 
in the SMC, the most metal-poor Cloud. 
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Table 1. Adopted logarithmic metal abundances in the SMC (on the scale in which loge[H]=12 and rounded to the nearest 0.05 dex). along with the source 
of the data. We consider the mean value (weighted by the number of lines) w hen abundances are g iven for different ions of the same element. The scaled solar 
mixture taking Fe as a reference is given for comparison (solar values from Asplund et al. 2009). The last row gives the resulting metallicity (assuming that 
the abundances of the elements not listed is 0.75 dex below the values of Asplund et al. 2009). 



Element SMC Scaled solar 

(Asplund - 0.75 dex) 



Source 



Reference 



C 

N 

O 

Ne 

Na 

Mg 

Al 

Si 

S 

Ar 

Ca 

Sc 

Ti 

V 

Cr 

Mn 

Fe 

Co 

Ni 



"zTxIcPT" 



7.50 
6.55 
8.00 
7.25 
5.60 
6.75 
5.60 
6.80 
6.40 
5.80 
5.55 
2.25 
4.30 
3.05 
4.90 
4.85 
6.75 
4.00 
5.30 



2.36 



7.68 
7.08 
7.94 
7.18 
5.49 
6.85 
5.70 
6.76 
6.37 
5.65 
5.59 
2.40 
4.20 
3.18 
4.89 
4.68 
6.75 
4.24 
5.47 



2.41 



B stars 

B stars/H II regions 

B stars 

H II regions 

Cool supergiants 

B stars 

B stars 

B stars 

B stars 

H II regions 

Cool supergiants 

Cool supergiants 

Cool supergiants 

Cool supergiants 

Cool supergiants 

Cool supergiants 

Cool supergiants 

Cool supergiants 

Cool supergiants 



1,2 

2, 3, 4, 5, 6, 7 

1,2,3,8 

4, 5, 6, 7, 9 

10, 11, 12, 13, 14 

1,2,3,8 

1 

1,2,3,8 

2 

4, 5, 6, 7, 9 

10, 11, 12, 13, 14, 15 

10, 11, 12, 13, 14, 15 

10, 11, 12, 13, 14, 15 

11, 12, 15 

10, 11, 12, 14, 15 

12, 15 

10, 11, 12, 14, 15, 16 

11, 12, 13 

10, 11, 12, 13, 15 



Key to references (for meaning of a cronyms see TablelAU — FlllKom et alj ( 
IVenneii & van der Hulslli20oA . r6l lKurt et aljh999h . r7l lGametllh999h . r8l1 
(1999), [12] L98 for MARCS models and spectroscopic gravities, [13] GW99, 



20001). [2] |Hunter et all <2005h. [3]|Hunter et alj J2007l) . [4] |Reveslh999l) . [5] 
JTrun dle et al. 1 20071). [9l lLebouteil ler et al. (2008). [10] Hilll h997h . [Ill Hill 
, [141 lVemitl999iu i5] RB89, ri6] .RQmaniello et aL 1200 A . 



As a first approach, we try to determine whether that source of 
discrepancy between theory and observation comes from the chem- 
ical mixture. For that purpose we have assembled for the SMC and 
the LMC a set of representative absolute abundance values for all 
elements significantly contributing to the opacity from a compre- 
hensive and critical examination of literature data. As the central 
bar is thought to have efficiently mixed the Interstellar Medium 
(ISM) through gas streaming, large-scale spatial chemical inhomo- 
geneities should not be a major issue in the case of young objects, as 
indee d confirmed by observations (e.g. H illl 19971 ; I Andrievskv et al.l 
uOOlh . The breakdown of these abundance data for all elements up 
to Z = 26 are shown in Fig.|AT]and|BT]for the SMC and the LMC, 
respectively. Some minor chemical species are not discussed be- 
cause of the lack of measurements, but their contribution to the 
global metallicity is negligible. 

The data are necessarily inhomogeneous and come from three 
independent sources: B stars, cool supergiants (note that they were 
chosen to cover a similar mass range as the B stars to ensure 
that they sample the same young population) and Hii regions. All 
these tracers suffer to different extents from various limitations and 
sources of errors (e.g. neglect of departures from LTE and evolu- 
tionary effects for the stellar objects, depletion onto dust grains and 
conversion from ionic to elemental abundances for the ionised gas). 
Although preference is given to the abundances of B stars, in a few 
cases we feel that data from other sources are more reliable and de- 
cide to use them instead (e.g. Ne data from Hii regions). For many 
elements, especially the heaviest ones, only data for cool stars are 
available. The NLTE corrections have been shown to be generally 
small and typically well wit hin ~ 0.2 dex, either for C, O and Na 
in LMC F-type supergiants JAndrievskv et al.ll200ll) or for O and 
Mg in SMC A-type supergiants jVennl 19991) . The LTE abundances 



of these metals in the evolved objects should therefore be relatively 
accurate, but this may not be true for other species (e.g. N in A- 
type stars: Venn 1999) . More generally, we give a higher weight to 
NLTE results and to the most recent studies based on higher qual- 
ity observational material and improved analysis techniques. We 
describe the available abundance data and the criteria used to select 
what we regard as the most probable present-day value for each el- 
ement in Appendices A and B. They are summarized in Tables [T] 
and|2]for the SMC and the LMC, respectively. 

As we aim to concentrate on the more metal-poor SMC, we 
note that the Fe abundance presented in Table [T] is in remarkable 
agreement with the average literature value for other young (< 1.5 
Gyr) SMC o bjects, that is < [Fe/H] > ^ -0.70 ±0.15 (gathered from 
ICarrera et a l. 2008; Parisi et al. 2009, 2010; De Propris et al. 201(3; 
Piatti et al. 2007a,b; I diart. Maciel & Costal2007i) . We also note the 
good agreemen t between our Z value and the mean value derived by 
iBuchleil ( 120081) . Indeed, we obtain Zsmc = 0.00236 and estimate the 
error to be ctsmc ~ 0.00079 (Zlmc = 0.00596 and ctlmc ~ 0.00232), 
while Buchler (2008) obtains Z = 0.0027 for the SMC (Z = 0.0046 
for the LMC). We determine the error on metallicity by a propaga- 
tion of the errors on individual element abundances. These errors 
are taken as the standard deviations of the different measurement 
distributions for each abundance (see Figures lATI and ISTi . When 
there are only a few data points for a given element, we use instead 
a fixed error of 0.15 dex, representative of typical errors given in 
the literature for abundances of MCs stars. 

Finally, we have determined the opacity profiles of a given 
stellar structure by considering, respectively, the AGS05 scaled so- 
lar mixture and the SMC mixture. The differences between the two 
opacity profiles do not exceed 3 %, while they are at most 6 % when 
the LMC mixture is considered instead. 
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Table 2. Adopted logarithmic metal abundances in the LMC (on the scale in which loge[H]=12 and rounded to the nearest 0.05 dex). along with the source 
of the data. We consider the mean value (weighted by the number of lines) w hen abundances are g iven for different ions of the same element. The scaled solar 
mixture taking Fe as a reference is given for comparison (solar values from Asplund et al. 2009). The last row gives the resulting metallicity (assuming that 
the abundances of the elements not listed is 0.35 dex below the values of Asplund et al. 2009). 



Element LMC Scaled solar 

(Asplund - 0.35 dex) 



Source 



Reference 



c 


8.00 


8.08 


B stars 


1,2,3 


N 


6.90 


7.48 


B stars 


4 


O 


8.35 


8.34 


B stars 


4 


Ne 


7.70 


7.58 


H 11 regions 


5, 6, 7, 8, 9 


Na 


5.90 


5.89 


Cool supergiants 


10 


Mg 


7.05 


7.25 


B stars 


4 


Al 


5.85 


6.10 


B stars 


1 


Si 


7.20 


7.16 


B stars 


4 


S 


7.10 


6.77 


Cool supergiants 


10, 11, 12 


At 


6.15 


6.05 


H 11 regions 


5, 6, 7, 8, 9 


Ca 


6.10 


5.99 


Cool supergiants 


10, 11 


Sc 


2.65 


2.80 


Cool supergiants 


10, 11, 12 


Ti 


4.70 


4.60 


Cool supergiants 


10, 11, 12 


V 


3.80 


3.58 


Cool supergiants 


11, 12 


Cr 


5.40 


5.29 


Cool supergiants 


10, 11, 12 


Mn 


4.95 


5.08 


Cool supergiants 


10, 11,12 


Fe 


7.15 


7.15 


B stars/Cool supergiants 


1,2, 10, 11, 12, 13 


Co 


4.65 


4.64 


Cool supergiants 


10, 11 


Ni 


5.85 


5.87 


Cool supergiants 


10, 11, 12 



"zTxKpy 



5.96 



6.05 



14 



Key to references — m lKometal] <200C|) . r2l lKom et al] <2002h . r3 llKomet alji2005h . f4l lHunter et al] i2008l) . r5l lRevesl {\99<k . r6l lLebouteiller et al.l 
<2008l) . [71 1 Verme il & van der Hulst (2002), [8] Peimbert ( 2003), r9l lGametll l ll999h . [101 AOl, [1 1] L98 for MARCS models and spectroscopic gravities, 
[12] RB89. [13l [Tmndle et aLt2007.) . ri41 .Romaniello et aH2Qoi) . 



3 STABILITY ANALYSIS OF MAIN SEQUENCE B-TYPE 
STARS FROM THE SMC 

We study the stability properties of MS B stars in the SMC. With 
the stellar evolutionary code CLES (Scuflaire et al. 2008), we have 
computed a grid of stellar evolution models with masses from 2.5 
to 18 Mo, a metallicity Z = 0.00236 and a He mass fraction fixed 
at y = 0.28. We have also computed some models with a lower 
abundance of He, Y = 0.26, which only affects the mass of models 
characterized by a given effective temperature (Teff) and luminos- 
ity (L). Convection is treated according to the mixing-length theory 
dBohm-Vitensdl 19580 with its parameter fixed at a^n = 1-8. The 
surface boundary conditions follow Eddington's law (r[T]) for a 
grey atmosphere. The convective core overshoot parameter is set to 
0.2 p ressure scale heights. The O PAL equation of state is selected 
from [Rogers & Navfonovl ( l2002h . Nuclear reaction rates are from 
NACRE (Angulo et al. 1999) with the revised "'N(p,7)'^0 cross 
section from Formicola et al. (2004). Concerning opacity, we have 
built new opacity tables for the specific SMC B star metal mix- 
ture (hereafter B-SMC mixture) w ith the last updated version of 
the OP codeu JBadnell et al.ll2005l) . Finally, we have performed a 
non- adiabatic stability analysis of these models with the code MAD 
(|Dupret2001). 

We do not find any SPB oryS Cep modes to be excited in mod- 
els at the SMC metallicity, Z = 0.00236. This is not surprising since 
opacity profiles with a scaled solar or a B-SMC mixture are very 
close and since results with solar mixtures already indicated no ex- 
citation at such a metallicity (Miglio et al. 2007a). In addition, we 
have made grids of models with the same physics, but with a metal- 
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Figure 1. Measurement of the instability of /3 Cep (left panel) low-order 
f = 0, 1 and 2 p modes, and SPB (right panel) high order ( = I and 2 
g-modes as a function of the period of the mode, log P (in days). Models 
are taken along the MS evolution phase for 12 and 4 M© with the B-SMC 
mixture and for different values of Z. 



licity of Z = 0.005 and Z = 0.01. At Z = 0.005, the B-SMC mixture 
only leads to excited SPB modes, while at Z = 0.01 we find a simi- 
lar instability domain as in the solar case for yS Cep stars and SPBs. 
Finally, we derive the minimum values of Z required to excite 



Pulsations in B stars of the Magellanic Clouds 



0.8 



0.6 



0.4 



0.2 




:::.^^JC-^^. 



Table 3. Set of parameters used to simulate an increase of opacity by means 
of Gaussians centred on tlie contributions of Cr, Mn, Fe and Ni. 



log T 



Figure 2. The solid line represents the OP normahzed opacity in a 12 M© 
model of a MS /3 Cep star. The individual contributions of Fe, Ni, Cr and 
Mn are represented and expressed in percentage of the total opacity. 



f} Cep and SPB modes in the SMC. We compute models of 4 and 
12 Mq, respectively representative of an SPB and /? Cep star by 
varying the metallicity from Z = 0.003 to 0.01 in steps of 0.001. In 
Fig.[T] we plot the normalized growth rate (;;): 



W 



JO \ dr \ 



(1) 



where R is the radius of the star. The represented 77 values come 
from the eigenmodes of models from the Zero Age Main Sequence 
(ZAMS) to the Terminal Age Main Sequence (TAMS). If 77 is pos- 
itive, the work (W) developed by the whole star is positive, mean- 
ing energy is efficiently transferred to drive the mode. It turns out 
that Z > 0.004 is required to excite SPB modes in the SMC, while 
Z > 0.007 is needed for/? Cep modes. 

In conclusion, the average B-SMC mixture and its correspond- 
ing metalUcity (Zsmc = 0.00236 + 0.00079) is marginally able to 
explain the presence of SPBs. But it is not able to explain the pres- 
ence of y6 Cep stars since it would imply that those stars have a 
metallicity significantly higher (at 6 cr) than the mean value for 
young SMC objects. 



4 CONSTRAINING THE IRON-GROUP OPACITY 

Already i n the eighties, an underestimation of the opacity was sug- 
gested bv lsimoni ( Il982h . He noticed that a revision of metal opaci- 
ties by a factor of 2 to 3 could explain the pulsations in;8 Cep s tars. 
New opacity computations by OPAL JRogers & Iglesiasll 19921 and 
references therein) confirmed it with a raise in the iron-group opac- 
ity by a factor of 3. More impro vements to stellar opaci ties con- 
tinued with the update of OPAL jiglesias & Rogerslll99q) and the 
work from OP JSeaton et al.ll994lBadnell et al.l2005h . Nowadavs. 



some discrepancies remain on the predicted iron-group bump be- 
tween OPAL and OP ( Seaton & Badnell 2004). Recent investiga- 
tions on monochromatic opacities of Fe and Ni reveal some dis- 
agreements between current opacity tables and new computations 
(e.g. Gilles et al. 201 1). Therefore, some room still remains for en- 
hancing opacities. Opacity calculations rely on the equation of state 



Element 
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5.07 
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Figure 3. OP opacity profiles in a 12 Mq /} Cep MS model (top panels) and 
a 4 Mq SPB MS model (bottom panels) with the SMC chemical mixture 
and metallicity. The non-modified opacity profiles are represented by black 
solid lines. The profiles increased by 10 and 50 % are drawn with long- 
dashed and short-dashed lines, respectively. They coincide with the normal 
opacity profile where no change is made, that is in cenfi'al layers (log T > 6 
) and towards surface (log T < 4.5). Opacity increases in the contributions 
of Fe, Ni, both Cr and Mn are presented from left to right. 



of stellar plasmas and on the detailed atomic transitions (lines) of 
elements that constitute the plasma. Provided the equation of state 
is correct, opacities can only be underestimated. Indeed, various 
approaches are used to implement transition lines in opacity codes, 
in particular with somewhat different restrictions on the number of 
lines. We consider this caveat as a gateway to possible underesti- 
mations. 

Hence we would like to constrain the opacity profile required 
to drive pulsations in MS B stars of the SMC. We proceed with an 
ad hoc local change in the opacity by considering different cases 
where the increase could be due to poor estimation of the opacity 
of respectively Fe, Ni, Cr and Mn. 

In Fig.[2l we show the normalized opacity profile (where k is 
the Rosseland mean opacity) in a 12 M© j8 Cep model and the con- 
tributions of each iron group element: in the opacity bump around 
log T = 5.3, Fe contributes to -65 %, Ni to -25% and Cr and Mn to 
less than 15%. Each contribution can be approximated by a Gaus- 
sian function centred at a value log T,, of width cTj (see Table |3]l- 
The changes in opacity are introduced by multiplying the original 
profile by a Gaussian with the same profile and an amplitude in- 
creased by A% (Table [3j. Let us note that a 50% rise in the Ni 
opacity corresponds to an increase of 50% of the total opacity at the 
temperature (log TiMi) at which Ni opacity has its maximum value. 
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Figure 4. Solid lines represent evolutionary tracks from ZAMS to TAMS of 
4 and 12 Mq models with a normal opacity. Long-dashed, dashed and short- 
dashed lines show evolutionary tracks of models with increases of 100% in 
the Fe, Ni, and Cr and Mn opacity, respectively. Filled circles indicate the 
model at X^ = 0.06 and 0.20 for the 4 and 12 Mq models, respectively. 



The Ni contribution alone is increased by more than a factor 2.5 
since its contribution to the total opacity in that region is ~ 25%. 

In Fig. [3] we illustrate, for given structures (log T , logp, X, Z) 
of respectively a /? Cep model (top panels) and an SPB model (bot- 
tom panels), the effects on the opacity profile of our ad hoc opacity 
modifications. Cr and Mn are treated together since their contribu- 
tions are small. As reported in Table |3] we consequently enhance 
the Mn opacity by one third of the increase of the Cr opacity since 
the contributions of Mn and Cr to the iron group bump follows ap- 
proximatively this ratio (see Fig.|2). Increases are always applied at 
the same temperatures, whatever the other quantities, on which the 
opacity explicitly depends, are. At a given temperature, the density 
(p) is not the same in an SPB as in a/3 Cep structure. Consequently, 
the maximum of the opacity increase in SPB models is shifted to 
slightly higher temperatures than inyS Cep models. 

With these modified opacities, we have computed new 4 and 
12 Mq models, representative of SPB and yS Cep stars, respectively, 
with Z = 0.00236 and physical inputs as detailed in Sect.[3] This in- 
duces changes in the density and in the temperature gradient (Vr), 
since a modification of the opacity directly affects VJ. As a higher 
opacity leads to expansion of the involved stellar layers, p decreases 
by 35-45% according to the element considered and the type of star, 
while vr is affected by changes of 10-15%. Yet, as illustrated in 
Fig.|4] effects in the Hertzsprung-Russell (HR) diagram are negli- 
gible since these changes take place only over a small fraction of 
the stellar radius. For a given evolutionary state, the stellar radius is 
increased by 1 to 2.5% either in/? Cep or SPB models. Whilst these 
changes do not affect the global properties of the star, they could 
however have some miscellaneous effect on mode excitation. 



4.1 Stability analysis 

We have checked the stability of these new models. The results are 
shown in Fig. [5] /? Cep modes are only excited if the opacity at the 
Ni maximum is increased by more than 50% (middle left panel). 
Changes at the temperatures of the Fe, and Cr and Mn maxima (top 
and bottom left panels, respectively) do not lead to /? Cep mode 



excitation, even for the 100% increase. Surprisingly, the more en- 
hanced the opacity in Cr and Mn is, the more stable the /? Cep 
modes are, as shown by the decrease in 77. SPB modes are excited 
for an increase > 30 % both in the Fe and Ni cases (top and middle 
right panels). We note that a 30% increase is very close to excit- 
ing SPB modes in the Ni case. No excitation of SPB modes is ob- 
tained in the Cr and Mn case, where the highest increases slightly 
strengthen the stability of SPB modes (bottom right panel). 

The T] parameter gives information on the global stability of 
a mode but it does not indicate which stellar layers contribute to 
the driving or damping of the mode. We thus introduce local quan- 
tities that describe more p recisely the k mec h anism, according to 
the formalism dev eloped in lUnno et al.l ( ll989l) . |Pamvatnvkhl ( 119991) 
and lPupretl 1 120011) . In particular, a stellar layer contributes to mode 
driving (resp. damping) if it produces a positive (resp. negative) 
work, i.e., -dW/d;?? > (resp. < 0). The latter quantity can be 
expressed as a function of the luminosity perturbation 5L: 



diy „ [ST d SL 

ex *K. < 

dm I T dm L 



(2) 



where 21 is the real part of the expression in brackets. We normalize 
W such that its value at the surface is equal to 3(a;), the imaginary 
part of the complex frequency of the mode. This equation is strictly 
valid only for radial modes, but remains a good approximation for 
non-radial modes of low angular degree, C. Indeed, corrective terms 
are generally <k d/dm{SL/L), due to the strong stratification in the 
radial direction. Assuming the diffusion approximation, we can fur- 
ther develop the right-side of Eq. (|2) as: 



T dm L 



Lr l^ST 6k dST/T 
dlnT 



5T' d [Lr I 5T 5k 

~ d^ [T \ T ~T 



Jr 



+A— + e(e+ 1) — )| 



,(3) 



where Lr is the radiative luminosity. The term 4^,/r -1- E{E + l)<f/,/r, 
where ^^ and f;, represent, respectively, the radial and horizontal 
displacements, is the geometric variation due to the oscillation. For 
radial modes, £i, = and t he term is thus sometimes called r (for 
radius) effect ( lBakeilll966h . It is one to several orders of magni- 
tude smaller than the other terms contributing to 6LIL and does not 
play any significant role. The term d(ST IT)ld\nT (hereafter the 
temperature gradient term) generally contributes to damping of the 
oscillations. Drivi ng due to the 4 6T IT term is sometimes called the 
7 mechanism (see lCox et al.ll 19661) , hence we will refer to it as the 
7 term. Finally, the term -Sk/k (hereafter the opacity term), with 
Sk being the opacity perturbation, is the largest of these terms and 
the main contributor to driving in the presence of an opacity bump, 
giving its name to the k mechanism. 

The presence of an opacity bump is not sufficient for the k 
mechanism to be efficient; there are two more necessary conditions. 
Firstly, the opacity bump must be located in stellar layers known as 
the transition region, i.e. where the thermal relaxation timescale, 
T|i,, is of the same order as the period of the modes, 11. Moreover, 
in that region, the eigenfunctions, in particular the pressure pertur- 
bation dP/P, must have a large enough amplitude and not vary too 
rapidly. Given that the value of T,h decreases monotonically from 
centre to surface in B stars, the periods of modes that are efficiently 
driven increase with the depth at which the opacity bump is located, 
which depends on T^g. This is why the frequency domain as well 
as the type of modes are clearly distinct between SPB and j8 Cep 
stars. Moreover, differences in the opacity profile for a given stellar 
model can thus alter the excitation of pulsations and the associated 
frequency range. 
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To analyse the effects of our ad hoc opacity modifications on 
the different terms of Eq.|3] we have chosen a 12 Mq model with a 
central hydrogen mass fraction X^ = 0.20 (as a representative fi Cep 
star), and a 4 Mq model before the second gravitational contraction, 
withXc = 0.06 (as a representative SPB). The location of these two 
models are indicated by filled circles in Fig. |4] In the next three 
subsections, we describe how the terms in Eq. |3](except for the r- 
term, which is negligible) react to modifications of opacity in the 
Fe, Ni, and Cr-Mn regions, respectively. 



4.1.1 The Fe case 

In Fig. [6] we show the variation inside the star of the work integral 
(panel a) for the fundamental radial mode of the p Cep model. In 
each panel, there are two different curves; dashed lines correspond 
to computations with a normal opacity, while solid ones correspond 



to the model with an opacity increased by 100% in the Fe domain. 
The derivative of the work integral (-dW/dra) is shown in panel b, 
and the different terms in the right-hand side of Eq. |3] are repre- 
sented in panels c to e. 

In region A, no work is performed since we are in the quasi- 
adiabatic part of the star where SS - 0, with SS being the entropy 
perturbation. On the other hand, region C has a very low heat ca- 
pacity and d/dm (6L/L) ^ 0. There, driving and damping contribu- 
tions balance to keep 6L/L constant. Finally in the stellar layers of 
region B, energetic processes that determine whether oscillations 
are efficiently driven or not, take place. 

It is known (see e.g. lUnno et alJI 19891) that the presence of 
an opacity bump due to an ionisation zone induces two tendencies 
in the derivative of the opacity. In the inner part of the ionisation 
zone, the derivative grows and contributes to driving the mode, as 
can be seen in region B (panel c) both in the normal and raised 
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Figure 6. From top to bottom are respectively represented the work integral 
(multiplied by 10*), the derivative of the work (i.e. the work developed by 
each stellar layer) as defined in Eq. 0, and the derivatives of first three 
terms of Eq. l|3]. The four latter quantities are multiplied by 10"* for a better 
visibility. The long dashed (resp. continuous) line represents the n= 1 , /" = 
mode of the 12 M© SMC /? Cep model computed with a normal opacity 
(resp. with a 100% opacity increase at the Fe peak) and X^ = 0.204. The 
grey shaded area represents the transition region. 



Table 4. High-order g modes excited in the SPB model with the Fe or Ni 
opacity increased by 50%. The minimal and maximal values of the period 
and of the ;; parameter found for these modes are indicated, rounded to 2 
and 4 decimal places, respectively. 



Element ( 



M) 



Fe 

Ni 



1.00-1.20 
1.63-2.47 
1.00-1.49 



0.0027-0.0221 
0.0012-0.0504 
0.0022 - 0.0539 



opacity cases. In the outer part of the ionisation zone, at lower tem- 
peratures, the opacity derivative decreases. This causes the opacity 
term to damp the oscillation in the region 4.7 < log T <5 (see panel 
c). As indicated by the negative value of W at the surface (panel a), 
the studied mode is not excited, either in the standard 12 Mq model 
or in the one with largely increased opacity. The reason can easily 
be understood: increasing the opacity around the Fe-maximum re- 
sults in a symmetric modification of the total opacity profile, and 
therefore driving and damping contributions are enhanced in a sim- 
ilar way. On the other hand, the moduli of the y and temperature 
gradient term contributions are enhanced in region B, leading to a 
stronger damping (panels d and e). As a consequence, the mode is 
stable even if the driving due to the k term increases. Actually no 
low-order p or g mode (of degree ^ = to 2) is excited even when 
the opacity around the Fe peak is increased by a factor of 2. 

In the SPB model with an opacity increased by 50%, 6 high- 
order (^ = 2) g modes are excited. Their ranges of periods and 
growth rates are given in Table |4] The contribution of the opac- 
ity term is sufficiently enhanced to drive these modes efficiently, 
whereas the contributions of the other terms do not show signif- 
icant differences between normal and increased opacities and do 
not influence the global stability of the mode. 



Figure 7. Same as Fig. |6] but now considering a model with an opacity 
increased by 100% at the Ni peak. 



4.1.2 The Ni case 

Nickel is of special interest. This element is typically 30 times 
less abundant than Fe (see Table [T]!, but its contribution to the 
iron-group opacity bump represents about half the Fe one. Since 
its contribution to the opacity peak is centred at higher tempera- 
tures, it could have a relevant role in the location of the blue border 
of the instability strip s as has already been demonstrated in sdBs 



or tne mstability stnp s 
JJeflFerv & Saiol2006l) . 



In thejS Cep model with a Ni opacity increased by 100%, we 
find the radial fundamental (n=l;^=0) and low-order (n=-l;^=l) 
and (n=-l;^=2) g modes to be excited. As the change is located at 
hotter temperatures, the driving and damping due to the k term are 
no longer affected in a symmetrical way, as was the case for Fe. 
In particular, raising the Ni opacity highly strengthens the driving 
in region B, as can be seen in panel c of Fig. [7] The moduli of the 
7 and temperature gradient term contributions are not amplified in 
regards to normal opacity conditions (region B, panels d and e). 
Therefore the effect of the ionisation zone on pulsations is biased 
towards mode driving when the Ni opacity is increased, and hence 
the/3 Cep fundamental mode is unstable (IV > in panel a). 

In the SPB model, increasing the opacity by 50% leads to 15 
dipole and 16 quadrupole excited high-order g modes, with peri- 
ods from 1 to 2.47 d (growth rates are indicated in Table |4]. The 
behaviour of an excited mode when the Ni opacity is increased is 
similar to the Fe case. The driving from the opacity term is in- 
creased so that the mode is excited while other terms are almost 
unaffected. 

However, while increasing the opacity by 50% at the Fe peak 
led only to 6 quadrupole excited modes, the same enhancement 
in the Ni region implies the excitation of 31 SPB modes (dipole 
and quadrupole ones). Furthermore, the instability domain is also 
affected. To illustrate this, we represent in Fig.|8]the high-order g 
modes found to be unstable during the MS evolution of the SPB 
model with the Fe and the Ni opacity raised by 50%, respectively. 
Firstly, the period domain is strongly increased between the Fe and 
Ni cases, in particular because ^=1 modes are never excited in the 
former case. Secondly, in the Fe case, unstable modes begin to be 
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Figure 8. Period as a function of the effective temperature for the C= 1 and 
1=2 excited high-order g modes in the 4 M© MS SPB models with the Ni 
and Fe opacities increased by 50% (represented in light/dark grey and black 
filled circles, respectively). The y axis is on a logarithmic scale. 



present at X^ = 0.15 (T^u ~ 15,900 K) while this hmit is at 
Xc = 0.55 (Teff ~ 18,000 K) in the Ni case, shifting to hotter 
temperatures the blue border of the instability strip of SPBs. 

As mentioned in Sect. [T] current theoretical models are not 
able to explain excited low frequency modes in Galac tic B-type 
pulsa tors, such as those found in 12 Lac or v Eri (Pamvatn vkh et al.l 
l2004h . The above described results motivate us to study the effect 
on the stability properties of a 10 Mq hybrid model with typical 
Galactic metallicity (Z = 0.014) of an enhancement of opacity by 
50% at the Ni peak. The periods of the unstable high-order g (SPB) 
modes are represented in Fig.|9]by black open (i=i) or filled (f=2) 
circles for the normal case and light {£=1} or dark (£=2) grey filled 
circles for the enhanced Ni case. The i=2 modes present a simi- 
lar instability domain. However the £=l domain is significantly in- 
creased in temperature and in period when the Ni opacity is raised. 
Hence there is a larger number of unstable modes with high period 
values, which are usually the ones not reproduced when modelling 
Galactic hybrids. Finally, the instability domain of/? Cep modes 
(which is not represented) remains unchanged from normal opacity 
conditions. 



Figure 9. Period as a function of the eft'ective temperature for excited 
high-order g modes in hybrid 10 M© MS models at Galactic metalhcity 
(Z=0.014). The light and dark grey filled circles represent the modes of de- 
gree t= 1 and f =2, respectively, for the stellar models computed with a Ni 
opacity increased by 50%. The black open and filled circles represent the 
modes of degree (=1 and (=2, respectively, for the stellar models computed 
with a normal opacity. The y axis is on a logarithmic scale. 
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4.1.3 The Cr and Mn case 

Increasing the Cr and Mn opacities does not allow either the p Cep 
or the SPB modes to be driven. More surprisingly, growth rates are 
weaker when the enhancements are larger (bottom panels of Fig. 

We can understand this result from the opacity term (panel c of 
Fig. not for the j8 Cep fundamental mode. The rise in opacity moves 
the driving component of the opacity term towards cooler tempera- 
tures outside the transition region of the fundamental mode. Since 
dLjL is constant in these low heat capacity regions, the driving is 
balanced by the y and temperature gradient terms. In such layers, 
every energy gain is directly lost as heat, inhibiting the conversion 
of the pulsation into mechanical energy. 

In the SPB model, with the increase of the opacity, the modu- 
lus of the opacity term contribution vanishes in layers where it was 
driving in the normal case. Hence the growth rate of the SPB modes 
is also decreased as the Cr and Mn opacities are enhanced. 



log T 

Figure 10. Same as Fig. |6] but now considering a model with an opacity 
increased by 100% at the Cr peak and 100/3 % at the Mn peak. 



5 DISCUSSION 

We have further investigated the consequences of increasing the Fe 
or Ni opacities in models where metallicities of Z = 0.003, 0.004 
and 0.005 are now considered. In the Fe case, no /? Cep modes are 
excited for Z = 0.003 even for an increase of 100%. However, p 
Cep modes are excited at Z = 0.004 and Z = 0.005 when 100% 
and 50% enhancements are applied, respectively. In the Ni case, 
for Z = 0.003, 0.004 and 0.005 we need increases of at least 100, 
100 and 50 %, respectively, to excite j8 Cep pulsations. Considering 
SPB models, we obtain excited high-order g modes at Z = 0.003 
for opacity increases of 20% and 30% in Ni and Fe, respectively. 
At Z = 0.004 an increase of 5% yields unstable modes both in the 
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Ni and Fe cases, as expected from results with the SMC mixture 
analysis in Sect.[3](see Fig.lTJ. 

We note that a determination of abundances for several Mag- 
ellanic B-type stars, including pulsator candidates, would allow us 
to assess if there is a clear deviation from the mean abundances 
-in particular the ones of iron group elements- we have adopted 
in this work. Though it would be valuable to establish the status 
of SPB and yS Cep candidates in the MCs, our work also presents 
interesting perspectives for Galactic hybrid B-type pulsators. Ad- 
ditional calculations for a Galactic model of a hybrid /? Cep/SPB 
star showed that, while the p Cep domain instability domain re- 
mains unchanged, the SPB domain is broadened when Ni opacity 
is increased. 



that current opacity data underestimate stellar opacity due to Ni by 
a factor <2. Therefore, we call for a reexamination of Ni opacity 
calculations. 
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6 CONCLUSION 

In this paper we search for an explanation for the presence of main 
sequence B-type pulsators in low-metallicity environments such as 
the Magellanic Clouds. Our investigation follows two approaches. 
First, we show that representative metal mixtures for B stars are 
not able to explain the presence of SPBs and especially of /3 Cephei 
stars in the more metal poor galaxy, i.e. in the Small Magellanic 
Cloud. 

In our second approach, we examine the effects of increasing 
the opacities of Fe, Ni, Cr and Mn in stellar models representative 
of the Small Magellanic Cloud. It is not possible to excite both 
SPB and /3 Cep modes with an enhancement of the opacity at the 
Fe peak, even for the largest increase. However, an increase of more 
than 50 % of the opacity profile at the Ni peak leads to the excitation 
of both SPB and p Cep modes. Finally, simultaneously raising the 
Cr and Mn opacities does not increase the probability of driving 
modes, instead it reinforces the stability of main sequence B stars. 

We also show that a modification of the opacity in the Ni re- 
gion moves the blue border of the SPBs instability strip to hot- 
ter temperatures and extends the frequency domain of pulsations, 
as expe cted from previous studies JPamvatnvkhlll999l : lMiglio et al.l 
l2007bh . 

Neither OPAL nor OP stellar opacities are able to reproduce 
the lowest frequencies (i.e. the highest periods) observed in Galac- 
tic hybrid B-type pul sators. A similar outco me in the modelling of 
sdB stars already led ljelferv & Said yOOg) to suggest an increase 
of the Ni contribution to opacity. Moreover, Ni (as well as Cr and 
Mn) atomic data used by OP to compute stellar opacities are not 
completely obtain ed ab initio. Instead , extrapolations from Fe re- 
sults are adopted JBadnell et al.ll2005h . These uncertainties in the 
iron-group opacities are at the origin of a multi-disciplinary collab- 
oration (see Gilles et al. 2011, for preliminary results), with the aim 
of comparing theoretical codes for opacity computations and carry- 
ing out experiments to measure opacities in conditions equivalent 
to those of stellar plasmas. 

Moreover, in the present study we find that an enhancement 
of the opacity in the region where Ni contributes could explain the 
presence of B-type pulsators in the SMC. Raising the Ni opacity 
also produces a larger period domain of excited modes in models 
of Galactic hybrid B-type pulsators, extending it in particular to 
higher periods (i.e. to lower frequencies) as required to explain ob- 
servational results. Already in the past, Simon (1982) made a plea 
for a reexamination of heavy element opacities. This was verified in 
the following decade when opacity due to the iron group increased 
by a factor of 3 for B star conditions. We hence conclude that B-type 
pulsators in the Magellanic Clouds and in the Milky Way suggest 
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APPENDIX A: SMC PRESENT-DAY ABUNDANCES 

We describe below the available abundance data in the SMC 
(LMC is treated in Appendix B) and the criteria used to select 
what we regard as the most probable present-day value for each 
element. They are summarised in Table [T] for the SMC. To ease 
readability, acronyms have been used in the following for the 
most cited references and are provided for convenience in Table 
All The absolute LTE abundances for B stars in NGC 330 of 



Acronym 


Reference 


AOl 


Andrievskv et al. (200r) 


GW99 


Gonzalez & Wallerstein (1999^ 


L98 


Luck etal. (199^ 


RB89 





Lennon. Dufton & CrowlevI ( 1200 3h are in poor agreement with 
other studies and are not considered further. 



Carbon Both iKom et alj ^2000') and ' Hunter et al] J2005l) sug- 
gest an NLTE value for B stars close to 7.5 dex. However, this 
estimate is only based on a handful of objects and must be cross 
checked against other sources. Despite the limitations affecting the 
carbon data from t he VLT/FLAMES survey (Hunter etal. 2007; 
iTrundle et al.l2007h . see Sect.|B] a close value is obtained from this 
much larger sample (~7.4 dex). On the other hand, all LTE studies 
of cool supergiants converge to a nearly identical abundance (Hill, 
Barbuy & Spite 1997; Hill 1999; GW9; RB89). This is perhaps 
fortuitous, as one may expect some depletion because of mixing 
(albeit negative, the NLTE corections should, however, be small; 
AOl). Equally pleasing is the good agreement with the data for H ii 
regions when corrected to first order for dust depletion (Fig. lAlb . 
suggesting that 7.5 dex is a robust estimate. 

Nitrogen Most B stars in the SMC show the products of CNO- 
cycle at their surface and are therefore N rich (as are a fortiori 
the cool s u pergia nts). This is the case fo r the stars studied by 
iKom et al .' (2000) and 'Trundle et al.l ( |200 7). Data for Hii regions 
cluster around 6.55 dex (Fig. lAlt . The consistency of these results 
and the negligible amount of nitrogen locked up in dust grains 
suggest that this estimate is reliable. Further more, it is identical t o 
the value found for the B-type star AV304 bv lHunter et al.l ( l2005h . 
Several B stars observed during the VLT/FLAMES survey also 
have upper limits consistent with such a low value dHunter et al.l 
l2007h . 

Oxygen All NLTE studies of B stars are in g ood agreement and 
the ipean value is adopted dKom et alj|200d ; [Hunter et alj|2005l 
|2007| ; iTrundle et all l2007h . With the exception of Hill (1999), 
our best estimate is close to the values found for A, F and K 
supergiants (Hill et al. 1997; GW99; Venn 1999 [NLTE]). This is 
also the case for the data for H ii regions (Fig. lAlt . 

Neon We base our abundance on observations of Hn regions 
dOarnetd [199^ iRevej IT99^; iKurt et al] 1 19991 ; iLebouteiller et all 
120081 ; IVermeii & van der Hulstil2002l) that only show a small scat- 
ter. Th e LTE neon abundance derived for AV304 by iHunter et al.l 
( I2OO5I) is obviously spurious (Fig. lAlt . 

Sodium There is no evidence for a strong Na enrichment arising 
from the dredge-up of Ne-Na cycled gas in SMC cool supergiants 
dnill 1997). We take the average of all the LTE values for yel- 
low/red supergiants (Hill 1997, 1999; Venn 1999; L98; GW99), 
but omit the anomalously high abundance reported for AVI 27 by 
RB89. 



Magnesium The numerous NLTE values in the literature for B 
stars indicate a value close to 6.75 dex without any discemable 
study- to-study differences dKom et alj l200d ; iHunter et all l2005l 
.20071 ; ITrundle etaLll2007l) . The (mostly LTE) data for evolved 
objects suggest a mean abundance ~0.1 dex higher, which may not 
be significant (Hill 1997, 1999; L98; GW99; RB89; Venn 1999 
[NLTE]). 

Alumi nium Only two B stars have NLTE abundances in the lit- 
erature dKom et al.ll2000l) , but the mean value is indistinguishable 
from the abundances found for F and K supergiants (Hill 1997, 
1999; GW99; L98 ). In contrast, the LTE abundance derived by 
iHunteret all d2005l) for AV304 is significantly lower. 

Silicon We take the average of the NLTE values found for hot star s 
dKom et alj|200d ; IHunter et alj|2005l , |2007| ; iTmndle et al.ll2007l) . 
but once again note a very good overall agreement with the data 
for cool supergiants (Hill 1997, 1999; GW99; L98; Venn 1999). 

Sulphur The adopted sulphur abundance is poorl y constrained 
and only relies on the NLTE value measured by IHunter et al.l 
(2005) in AV304. L98 obtained a LTE value -0.16 dex higher in 
one Cepheid, which is consistent within the errors. The (slightly) 
lower value inferred from H 11 regions (Fig. lAlt also provides some 
support for our estimate, as some dust depletion is expected. 

Argon All studies of the ionised gas in Hii regions yi eld very 
simil ar results and we assume the mean value here dGametn 
,1999'; 'Reyes' '1999"; 'Kurt et alJ 1 19991 ; ILebouteiller et alj I2OO8I ; 
Vermeil & van der Hulst 2002). 

Calcium, Scandium and Titanium We take the median of the 
abundances derived for cool stars (Hill 1997, 1999; GW99; L98; 
RB89; Venn 1999). 

Vanadium There is a large spread in the published values 
(Fig.lAlJ. We average the values of Hill (1999), L98 and RB89. 
The abundance of this element is uncertain. 

Chromium We consider the typical value found for cool super- 
giants (Hill 1997, 1999; L98; RB89; Venn 1999). 

Manganese The values found by L98 and RB89 have been used. 

Iron Compared to the LMC, much fewer B-type stars have had 
their iron content determined and the scatter from one study to 
another is also slightly larger: Hunter et al. (2005; 6.63 dex in 
LTE), Kom et al. (2000; 6.82 dex in NLTE) and Tmndle et al. 
(2007; 6.93 dex in LTE). We prefer to rely on the numerous 
observations of yellow/red supergiants that suggest an abundance 
-6.75 dex (Hill 1997, 1999; L98; RB89; V enn 1999; R omaniello 
et al. 2008), but note that the NLTE value of lKom et al] dJOOOt) for 
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two B stars is only 0.07 dex away from this value. 

Cobalt We use tlie values of Hill (1999), L98 and GW99, but omit 
the star with an anomalously low abundance in the sample of Hill 
(1999). 



from convective mixing in LMC cool supergiants (AOl). We 
adopt the mean NLTE abundance of AOl for F supergiants, which 
compares favourably with the values of L98 for Cepheids if NLTE 
corrections of the order of -0.1 dex are applied (AOl). The value 
for one star of RB89 is discrepant. 



Nickel The values of Hill (1997, 1999), GW99, L98 and RB89 
have been used. 



APPENDIX B: LMC PRESENT-DAY ABUNDANCES 

As for the SMC, we present here a compilation of abundances for 
the LMC and the criteria used to select what are the most probable 
present-day value for each element. They are summarised in Table 
[2l Acronyms of Table lAll are still used in this section. 

Carbon We adopt the NLTE abundances derived by Korn et 
al. (2000, 2002, 2005). Carbon abundances are available for the 
large number of stars obse rved in the course o f the VLT/FLAMES 
survey of massive stars ( IHunter et al.l l2008h . However, because 
these results are mainly based on lines of which reliability is 
doubtful (e.g. Cii /14267) w e do not consider these results here . 
The few LTE abundances of iRolleston. Trundle & DuftonI ( l2002n 
are in poor agreement with other studies. The slightly lower 
abundances found for cool supergiants (~0. 15 dex; RB89; AOl 
[NLTE]), if physical, might be explained by evolutionary effects. 
The values for Hii regions are in very good agreement (Fig. IBlt . 
Only a small amount of carbon is believed to be trapped in dust 
grains (~0. 1 dex in the Orion nebula according to Esteban et al. 
1998). 



Magnesium We c hoose the mean value of 7.05 dex found by 
iHunter et alj ( 120081) . Other studies based on smaller samples of 
early-type stars give very s imilar results (IKom et al.lbOOOl . |2002| . 
I2OO5I ; iRoUeston et al j[20o3) . Mean LTE abundances between 6.88 
and 7.44 dex have been reported for cool supergiants (AOl; L98; 
RB89). 



Aluminium The LTE value of lRoUeston et all j 19961) for PS34-16 
is on the low sid e (Fig. IBlt . We adopt the NLTE abundances of 
iKom et alj <2000l) . Higher values are suggested from LTE studies 
of F supergiants (-6.1 dex; AOl; L98). 

Silicon The Si abundances of a large sample of B stars have 
been derived in t he framework of the VLT/FLAMES survey 
JHunter et alj200 6) and we use these values her e. This mean valu e 
is in agreement with the LTE estimates of Roll eston et al.l ( l2002h . 
while it lies in between the results for other B stars in NLTE 
JKom et alJl2"oO0L I200I I200I) and red supergiants in LTE (AOl; 
L98; RB89). The abundances from Hn regions are unreliable 
because of severe dust depletion. 



Sulphur The LTE abundance estimated bv IRolleston et al.l 12003) 
clearly appears underestimated (Fig. IB 11 . We adopt the mean 
of the few LTE values reported for evolved objects (AOl; L98; 
RB89). The significantly lower abun dance found from H n reg ions 
likely arises from dust depletion (e.g. lLebouteiller et aLl2008h . 



Nitrogen Nitrogen abundances show a very wide spread in B stars 
in the LMC, with most of them exhibiting a strong excess as a 
result of mixing. We adopt the baseline abundances of lHunter et al.l 
( |2008[) coiTesponding to the lowest values found in their sample. 
This is consistent with the v alues of Korn et al. ( 2000, 2002, 
2005). The stars analysed by IRolleston et al.l ( I2OO2I) are clearly 
enriched. The data for Hii regions show a lar ge scatter (Fig. IB 11 . 
but our estimate is in excellent a greement with lOamettl d 1 9991) and 
IVermeii & van der Hulsg 00021) . Nitrogen is not expected to be 
significantly depleted in the gas phase. 



Oxygen The oxygen abundance is taken from lHunter et al.l ( l2008h 
and is in very good agreement with the values of Rolleston et 
al. (2002; LTE) and Korn et al. (2000, 2002, 2005; NLTE). This 
value is reasonably consistent with the NLTE values of AOl for F 
supergiants (~8.5 dex) and the results for Hii regions (Fig. IB It . 
Only a small amount of oxygen is believed to be t rapped in dust 
grain s (~0.08 dex in the Orion nebula according to lEsteban et al.l 
Il998l) . 

Neon We base our abunda n ce on observation s of H n regions 



JGarnetll IT999I; iRevesI [l999l; jLebouteiller et all I2OO8I; IPeimberll 
l2003l ; IVermeii & van der Hulsll2002l) . There is a good concordance 
between these values a nd, as a n oble gas, neon is not embedded 
in dust grains. Rollest on et al.l 11996) derived the LTE neon 
abundance in LH 104-24, but this value is discrepant and very 
likely grossly overestimated (Fig. IB It . 

Sodium There is no evidence for a strong Na enrichment arising 



Argon No stellar studies exist and we base our results on ob- 
serva ti ons of Hn reg i ons, which are in c lose a g reement dGametd 
I999I; iRevej 1 19 99V Leb outeiller et al.l I2OO8I ; IPeimberd I2603F 
Vermeii & van der Hulst.20021) ^ 



Calcium We only use the values of AOl and L98, as the values 
of RB89 based on lower quality data seem slightly discrepant 

(Fig.ED. 

Scandium There is a large spread in the results of the different 
studies (AOl; L98; RB89), but we take the mean as a representative 
value. 

Titanium Excellent agreement is found between the various 
studies (Fig. lBlt and we adopt the mean value (AOl; L98; RB89). 

Vanadium We average the values of L98 and RB89 and omit the 
discrepant value found for one star by AOl (Fig. lBlt . 

Chromium and Manganese AOl find on average slightly lower 
values than L98 and RB89 for both elements. As the distributions 
overlap, however, we simply take the mean value. 



Iron We ignore the LTE value of IRolleston et alj ( l2002h for one 
star based on 2 lines differing by up to 0.6 dex. The mean values 
for the other B star studies are in fairly good agreement: Korn et al. 
(2000; 7.09 dex in NLTE), Trundle et al. (2007; 7.21 dex in LTE) 
and Korn et al. (2002; 7.3 2 dex in LTE). Con sidering the smallness 
of the NLTE coiTections JKom et al.ll2000l) . we can compute the 
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Figure Al. SMC abund ances. B stars black: Hunter et al (2005, 2007) and'Xrundle et al.' feOOVl), rerf: lKomet"ZI <2000h Cool supergiants magenta: L98 and 
iRomaniello et alj i2008l) . blue: GW99 cyqwilVennl h999l) update d bylv enn & Przybilla (2003) fo rN, green: RB89,ye//ow : Hill et al. (1997) and Hill (1997, 
1999) H II regions dotted line: Gamett fl99y), short-dashed /mei lRevea 1 19991) . long-dashed /mei lLebouteiller et alJJ200 }|). dot-short dashed line: Kurt et al 
(1992), dot-long dashed /;'«e: .Vermeil & van der Hulst t2002h . The vertical, red ticks at the top of each panel indicate the solar abundances of .Asplund et al, 
(2009^ - 0.75 dex. 



grand average of all these values to obtain ~7.21 dex. On the other 
hand, the LTE analyses of the F-type supergiants suggest ~7.14 
dex (AOl; L98; RB89; Romaniello et al. 08). We assume 7.15 dex 
as a compromise. Iron is dramatically depleted in dust grains, such 
that Hii regions do not offer any constraints (Fig. IB It . 

Cobalt We take the mean of the values of AOl and L98. 

Nickel We take the mean of the values of AOl, L98 and RB89. 



Pulsations in B stars of the Magellanic Clouds 15 



16 S. Salmon et al. 




10 



5 
4 
3 
2 
1 


5 
4 
3 
2 
1 


5 
4 
3 
2 
1 


5 
4 
3 
2 
1 




7 


7.5 




8 


8.5 


■■■■!■ 

"Ne 


■■■I| 


1 


'!' 1 


■■■!■■ 


, . . . 1 , 


. . . 1 


1 

; 


.!. 1 


,,I,,T 



7 7.5 


8 


8.5 


1 1 1 1 1 1 1 1 
- Al 


. , 1 . . , 


1 1 1 1 1 


1., ..il, 


III 


1 1 1 1 


, , ..1 



5 5.5 6 


6.5 7 


1 1 1 1 1 1 1 1 1 i| 1 
-Ca 

"..I..II1 Il 


1 1 1 1 1 1 1 1 1 

,,iii,,,r 



5 5.5 


6 6.5 7 


- V 


: "Ill 


1.11 ji...,..: 




6.5 7 7.5 8 
log £ [dex] 



^^TTTT 




16 - 



I"" I 



W....I 



6.5 7 7.5 8 8.5 



" ' I " T I " " I " " I 




I ■■■■■ 



20 



5.5 6 6.5 7 



5 E- 




T 



I T ' I " " I "d 



15 j- Si 
10 I- 



■ . I ■■^. >.iil 



Ml 



I'liL i-ii I I I I ■- 



5 
4 
3 


6.5 7 


■ 7.5 8 


1 1 1 1 1 1 1 1 1 . 1 .. 1 . 1 . . 1 . 1 . 
"Sc 


2 
1 

n 


Liiil 


.l.li.i.,.,: 




4.5 5 5.5 6 6.5 




3.5 4 4.5 5 5.5 
log e [dex] 



12 



7.5 8 8.5 9 



-I 10 - 

8 - 

6 - 

4 - 

"1 2 - 







Mg 



'|"l ""l_L 

A 



Jbli 



6 6.5 7 7.5 8 




3.5 4 4.5 5 5.5 




4 4.5 


5 5.5 


6 


1' 

"Ni 

1 


1 1 1 1 1 1 1 1 


' 1 


1.., . 


iril 


1.: 


.1 



5 5.5 6 6.5 7 
log c [dex] 



Figure Bl. LMC abundances. B stars black: iHunter et al.l i200d) and[Trandle et al.' ('2007"), /j)/Me:'Rollesto n et alj (20o3 ). red: Korn et al. (2000, 2002, 2005) 
Cool superg iants magenta: L9 8, cyan: API, green: RB 89, ye/tow: [Romaniello et al. (2008) H II regions dotted line: Gamett (1999), short-dashed line: 
iRevea 419991) . long-dashed /me: lLebouteiller et al.l 420081) . dot-short dashed line: Peimbert (2003), dot-long dashed /mei lVermeii & van der HulslU2002h . The 
vertical, red ticks at the top of each panel indicate the solar abundances of iAsplund et al.i l2009i) - 0.35 dex. 



